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Abstract 



We review the mechanism of weak decay of hypernuclei, with em- 
phasis on the non-mesonic decay channels. Various theoretical ap- 
proaches are discussed and the results are compared with the available 
experimental data. 



PACS: 21.80.-Fa; 13.75.Ev; 25.40.-h 

1 Introduction 

The main decay modes of A-hypernuclei are the so-called mesonic and non- 
mesonic decays (for a recent review on the argument see ref. jT|): the former 
also occurs for free A-hyperons: 



The experimental ratio of the relevant widths, T^^'^/r^o'^ ~ 1.78, together 
with the measurements of A polarization observables lead to the formulation 
of the A/ = 1/2 rule on the isospin change of the system [from simple 
Clebsch-Gordon coefficient analysis, the latter would predict F^f^'/r^o*^ = 2]. 
This rule is based on experimental observations, but its dynamical origin is 

*Presented at XXVIII Mazurian Lakes Conference, Krzyze, Poland, August 31- 
September 7, 2003 




(pfree/pfree ^ o.639) 
jpfree/pfree ^ g.SSS) 



(1) 



not yet understood on theoretical grounds. The Standard Model does not 
support it and many non-perturbative effects could be responsible for the 
measured enhancement of the AI = 1/2 transition amplitude. 

The mesonic decay can also take place in A-hypernuclei, however the 
Pauli principle tends to disfavour the produced final state, with an emitted 
nucleon of momentum ~ 100 MeV/c (from a A at rest), well below the 
maximum level of occupied states in the nucleus. Of course this argument 
is only a qualitative one: indeed the mesonic decay is an important fraction 
of the total decay rate in light-medium A-hypernuclei. This outcome stems 
from several facts: i) the hyperon momentum distribution in the nucleus, ii) 
the medium attraction on the emitted pion, which lowers its mass at fixed 
momentum and distorts the pionic outgoing wave, in) the strong reduction 
of the local Fermi momentum at the nuclear surface. [2] 

Notably the mesonic decay can be used to extract information on the 
pion~nucleus optical potential (in a complementary way with respect to 
pionic atoms and low energy tt- nucleus scattering experiments). 

The A-hypernuclei non-mesonic decay occurs via the weak interaction of 
the A with one nucleon: 

AA^ NN (one-body induced decay) 

or with a pair of correlated nucleons: 

KNN NNN (two-body induced decay). 

These processes are mediated by the exchange of a meson (including the 
K, K* strange mesons), as illustrated in Fig. ^ 
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Figure 1: One-nucleon and two-nucleon induced A decay in nuclei 
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The decay width associated to the various non-mesonic processes can be 
denoted as follows (n indicates a neutron, p a proton): 





An 


— > nn 


(2) 




Ap 


— > np 


(3) 


r2 : 


ANN - 


NNN 


(4) 



whence the total non-mesonic width is Fnm = 1^1 + ^2 = + + The 
total hypernuclear decay width is then: 

Tt = Tm + Tnm (5) 

where Fm = T^- + F^o is the mesonic width. 

The non-mesonic mode can occur only in nuclei, and therefore it is the 
only source of information on the AA^ — > NN weak interaction. 

Concerning the kinematics of the process, it is worth noticing that if 
one assumes the A at rest (Q-value of about 176 MeV~ m\ — nrij^i) and 
an equal share of the available energy between the two (respectively, three) 
final nucleons, the momenta of the emitted nucleons are expected of the 
order of pN — 420 MeV/c for the one-body induced decay (respectively, 
Pn — 340 MeV/c for the two-body induced decay). Hence we do not expect 
a significative influence of the Pauli principle on these decay channels. 

An interesting feature of the experimental decay widths of A-hypernuclei 
is the approximate stability of this quantity in passing from medium-light 
to very heavy systems: the mesonic (Fm) and non-mesonic (Fnm) partial 
widths depend upon the mass number in such a way as to compensate each 
other. This saturation property is clearly related to the short range of the 
AN weak interaction. [3 1^ 

2 Theoretical models for A- hypernuclear decay 

The starting point for all methods described below is the weak effective 
hamiltonian for A — > vrA^ decay: 

'HZN = iGml'^^{A + B^^)f-$^^K, (6) 

the parameters of which are fixed on the free A decay: G = 2.211 • 10^^/m^, 
A = 1.06 (PV amplitude) and B = —7.10 (PC amplitude). To enforce the 
A/ =1/2 rule the hyperon is assumed to be an isospin spurion with I = 1/2, 
h = -1/2. 
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In the non-relativistic limit, one can then express the free A decay width 
as follows: 




where Ca = 1(2) for F^o (F^-, respectively) and the decay occurs both 
through parity violating (the s-wave amplitude, S = A) and parity conserv- 
ing (the p-wave amplitude, P = uItj^B /{2mN)) terms. 



2.1 The wave function method 

An expression similar to for the mesonic width of hypernuclei can be 
obtained by explicitly taking into account the hyperon (i^a), pion (cj),^) and 
nucleon {(Pn) wave functions inside the nucleus [210 El : 



(8) 



X 5^ 



dr(t)K {f)(f)Tr {q, r)(f)*^ (f ) 



p2 

ml 



df(j)A (f ) V(/>7r (g, r) 4>*j^ {r] 



The pion wave function corresponds to an outgoing wave solution of the 
Klein-Gordon equation in the presence of a suitable vr-nucleus optical po- 
tential, K)pt, while the A wave function can be derived within a variety of 
hypernuclear shell models, whose parameters are typically determined by 
comparison with the available spectroscopic data. 

Turning now to the evaluation of the non-mesonic width, one needs 
an explitic model for the AA^ — > NN weak transition: the latter is usu- 
ally described in terms of the exchange of virtual mesons belonging to the 
pseudoscalar (vr, 77 and K) or to the vector (p,w and K*) octets. The most 
important (from an heuristic point of view) component of the AA — > NN 
transition potential is associated to the exchange of a pion: 



V^{q) = -Gm. 



2 gNN-n 
2mN 



~^ 2m 



0"2 



(p + mj, 



T2 



(9) 



with fh = {m\ + TTiTv) /2. 

This is the only component of the AA — NN meson exchange poten- 
tial, for which the couplings of both the weak and strong vertices are well 
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constrained by the existing phenomenology and experimental data. The 
exchange of heavier mesons, which cannot be produced in the free- A weak 
decay, is subject to theoretical uncertainties and to a somewhat large in- 
determination of the model parameters (typically couplings and masses), 
in spite of the existence of various theoretical schemes which, in princi- 
ple, should allow to fix them on a firm basis. Nevertheless the influence of 
AA'^ NN processes mediated by heavier mesons on the calculated non- 
mesonic decay widths appears to be important in order to reproduce the 
available experimental data. 

A sort of "minimal model" , which has been often employed in the liter- 
ature .7^, accounts for the exchange of pions and p-mesons, together with 
phenomenological, g-dependent short range correlations, usually parameter- 
ized in the Landau-Migdal form (see, for details, ref. ^). 

In the framework of the wavefunction method, the one-body induced 
non-mesonic decay width takes the form: 

^'=IWfI (§3 2-^(E.C.)El-^(Pl,P2)l^ (10) 

where 

M{pi,P2) = {^R;N{pi)N{p2)\fAN-.NN\'^H) (H) 

is the transition matrix element between the initial hypernuclear state \^h) 
and the final state with two outgoing nucleons, while 5(E.C.) guarantees en- 
ergy conservation. Besides the above mentioned uncertainties in the OME 
(One-Meson-Exchange) transition potential, one faces here also the compli- 
cations associated with a strongly interacting many-body system. 

We will not consider here further details of this method, however it is 
worth recalling that the wave function approach is considered to be the most 
appropriate one for the calculation of light A-hypernuclei mesonic widths; 
it is, instead, less accurate for the non-mesonic widths and certainly less 
practical for the heavier systems. 

2.2 The polarization propagator method 

This method obtains the hypernuclear decay width through the evaluation 
of the A self-energy inside the nuclear medium |51 13 13 El • 

rA = -2ImSA. (12) 

The starting point is again the weak effective Hamiltonian which 
provides the relevant vertices already in the simplest process illustrated in 
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Figure 2: Lowest order terms for the A self-energy in nuclear matter. 



diagram (a) of Fig. [2 which corresponds to the free A decay. More generally, 
a few contributions to the A self-energy in the medium are illustrated in 
diagrams (b)-(h) of Fig.|2 where the pion propagator is typically dressed by 
the strong interaction with particle-hole (ph), A-hole (Ah) and two particles- 
two holes (2p2h) intermediate states. Further contributions arise from the 
medium modifications on the nucleon propagator, which has origin in the 
AvrA^ weak vertex. 

Formally the A self-energy can be written as: 



Sa(A:) = 3i{Gmlf 



4 J 52 + ] Fl{q)GN{k - q)G^{q) , (13) 



where k is the A momentum inside the nucleus, Gi\f{p) and GT^{q) are the 
nucleon and pion propagators (in nuclear matter): 

-A;^) , e{kF-\v\) 



Gn{p) 



P 



+ 



and 



po - En{p) -Vn + ie po- En{p) - Vn 



ql - g2 



mt 



Kid) ■ 



(14) 
(15) 
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In equation (|14() the medium effects on tlie nucleon propagation are plie- 
nomenologically embodied in the nuclear binding potential Vat, while in ()15|) 
the analogous effects on the pion appear in the pion self-energy, illustrated 
above. 

More explicitly, by inserting ()13|1 into p2|l. one can express the A- 
hypernuclear decay width as: 

TA{k, p) = -%{Gmlf I jPyreH k-q\ -kp) 



where: 



. (2^)3 

xe(/co - EN{k -q)- yjv)Im [a{q)\^^^^_j^^^^_^_y^ (16) 



I Pl{q)UL{q) , . P^{q)UT{q) 

^l-Vz.(<?)f/L('Z) 1 - 1/T(g)f/T(<z) ^ ^ 

In the above, to perform a realistic calculation, the effective interactions 
5, Pl, Pt,Vl, Vt include vr- and p-exchange plus short range repulsive cor- 
relations. Specifically: Vl, Vt are the (strong) p-h interaction and include a 
Landau parameter g'{q), which embodies the NN short range repulsion. S, 
Pl and Pt correspond to the lines connecting the weak and strong hadronic 
vertices and contain another Landau parameter, g'/^{q) [a priori different 
from g'], which is intended to parameterize the strong AA^ short range cor- 
relations. 

Besides the interaction lines, the other key- ingredient of l|17|) are the 
longitudinal and transverse (with respect to the pion momentum q) polar- 
ization propagators, Ul and Ut- They contain the Lindhard functions for 
p-h and A-h excitations and the irreducible 2p-2h polarization propagator: 

UL,Tiq) = UP\q) + U^\q) + Ulf\q) (18) 

The imaginary part of a{q), needed to evaluate La, develops various 
contributions, namely: 

UL,T{q) _ Im UP\q) + Im U^%) + Im Ulf\q) 

""^-VlADUlAD \i-Vl,t{<i)UlAq)? ' 

Indeed the three terms represent different decay mechanisms of the hyper- 
nucleus: 

Tm oc Im U^^ (part of mesonic width) 
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Fi cx Im U^^ (non-mesonic, one-body induced decay width) 

r2 oc Im u'^P'^^ (non-mesonic, two-body induced decay width 

and additional part of mesonic width) 

While the p-h and A-h polarization propagators are well known and can 
be analytically evaluated il2j, the 2p-2h polarization propagator, even in 
the non-relativistic limit considered here, demands quite a computing effort 
to be fully determined. On the other hand, according to the above relations, 
it is needed for estimating the two-body induced decay width. 

In this context two approaches have been utilized till now for the evalu- 
ation of Uf^^^^{q): 

A. Phenomenological model. 

In this case one takes into account mainly the phase space available 
to 2p-2h excitations in nuclear matter and determines the relevant 
entity of the 2p-2h propagator through its connection with the phe- 
nomenological TT-nucleus optical potential. The relation with specific 
hypernuclei is then obtained by implementing the local density ap- 
proximation. 

B. Microscopic calculation. 

This approach proceeds through the path-integral formulation and 
provides the 2p-2h propagator within the so-called One Boson Loop 
(OBL) approximation, which embodies a rich variety of perturbative 
contributions to Ul^^^{q). 

In the following we will consider only the phenomenological approach and 
we refer the reader to refs. ^2 Q ^oi the details and the results of the 
microscopic approach. 



2.3 The phenomenological 2p2h propagator 



In the (goj Q) region where the p-h and A-h excitations are off-shell, the fol- 
lowing rela 
Vopt holds: 



lowing relation between jj"^"^^ and the p-wave pion-nucleus optical potential 



j2p2hi 


l-^5L(g)f/L(g) 



2(7oKpt(g) . (20) 



At pion threshold the latter can be parameterized as: 

2goKpt(go ^m^A^ 0;p) = -47rg V^<^o (21) 
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where Cq is a complex number extracted from the experimental data on 
pionic atoms its "proper part" (using g' = gii^) = 0.615) turns out to 
be: 

Co* = (0.105 + i0.096)/m^ . 

Hence one obtains the following parameterization of the proper 2p-2h po- 
larization propagator in the spin-longitudinal channel at pion threshold: 



r^I^F^Ul^'^qo ^ m^,q^ 6; p) = -47rg V'Q* ■ (22) 



Further, to obtain the general dependence of u"^^ upon {qo,q), one consid- 
ers the phase space available for the real 2p-2h excitations: 

Piqo, q; p) oc I ^ Im U^'^ (| + k; p) 

xlmU^^ (I - (f + A:o) 9 (f - fco) , (23) 

which amounts to neglect the energy and momentum dependence of the p-h 
interaction. Finally the imaginary part of u"^^ will be written as: 

Im C/f^^go, P) = f(g°'|^\ lm C/ff (m., 0; pefr) (24) 

with PefT = 0.75/9. 

The above formulas refer to homogeneous nuclear matter and provide, 
through equations (jlOj) and H17|) , the hypernuclear decay width for a A with 
momenum k embedded in a constant nuclear density p. The corresponding 
decay width in finite hypernuclei can be obtained by applying the local 
density approximation (LDA). 

The latter amounts to consider a local Fermi momentum 

fcF(r) = |-vrV(r)| , (25) 
which is defined in the Thomas- Fermi approximation, as follows: 

eF(r-) + VN{f^ ^ ^ + VN{r) = 0. (26) 
2m7v 

The decay width in finite nuclei is then obtained by: 

rA= / d^|^A(^)|2rA(fc) (27) 
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where: 




(28) 



3 Theory versus experiment 

In the phenomenological approach for the polarization propagator we have 
employed the customary Fermi distribution of the nuclear density: 



with R{A) = 1.12^1/3 - 0.86y4-^/3 fm and a = 0.52 fm. 

The A wave function is obtained from a Woods-Saxon A-nucleus poten- 
tial, which exactly reproduces the first two single particle eigenvalues {s and 
p A levels) of the hypernucleus under analysis ,10,. 

Short range correlations {g' , g'^^) are fixed to get agreement with exper- 
imental data: indeed the Landau parameter g' has been widely used in the 
literature, in connection with spin-isospin nuclear responses as measured in 
charge exchange reactions and inelastic electron scattering. In this context 
the customary value, used within the RPA framework and appropriate to 
fit the data is g' = 0.6 H- 0.7. 

The available data from A-hypernuclei decay seem to require a somewhat 
larger value of g' (0.8), together with g'j^ ~ 0.4. This is not necessarily in 
contrast with previous phenomenology, since the RPA-type correlations are 
now applied to a richer polarization propagator, which also includes the 
2p-2h excitations. 

In Fig. 121 we illustrate the results obtained ^U] with the polarization 
propagator method (together with LDA) for the decay widths of various 
A-hypernuclei. In addition to the total decay rates, the mesonic and non- 
mesonic partial rates are shown, the latter being separated into the one-body 
induced (Fi) and two-body induced (F2) decay rates. The total decay rate 
appears to be rather constant with the mass number, at least for A > W 
(saturation property), as a result of a compensation between the mesonic 
decay and the non-mesonic one. 

Altogether the results of the theoretical calculation appear to be in good 
agreement with the available experimental data, both on the total rates and, 
when available, with the non-mesonic rates. Similar outcomes, though not 




(29) 
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Figure 3: Partial A decay widths in finite nuclei as a function of the nuclear 
mass number A. The experimental data are taken from Refs. jl4l 1151 116j 
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reported here, were obtained with the microscopic calculation of the 2p- 
2h polarization propagator, a fact which demonstrates how the theoretical 
description of these processes is well founded. 



4 The Tn/Tp puzzle 

The main problem concerning the weak decay rates is to reproduce the ex- 
perimental value for the ratio Tn/Tp between the neutron- and the proton- 
induced widths: 

An nn and Ap — s- np (30) 

Theoretical calculations underestimate the central data for all considered 
hypernuclei: 



< 2 (31) 



One should keep in mind that, up to now, the data on the separate 
rates are limited and not precise enough, due to the difficulty in detecting 
the products of the non-mesonic decays, especially neutrons. The present 
experimental energy resolution does not allow to identify the final state of 
the residual nuclei in —>■ ^~^Z + nn and —>■ ^~^(Z — 1) + np. 

In the OPE approximation, by assuming the AI = 1/2 rule in the A 
7r~p and A — > ir^n free couplings, many different calculations give small 
ratios, in the range 0.05 ^ 0.20. 




OPE 

~ 0.05 H- 0.20 (32) 



for all the considered systems. 

For pure AI = 3/2 transitions the OPE ratio can increase up to about 
0.5. However, this assuption would be inconsistent with the fact that the 
OPE model with AI =1/2 couplings well reproduces the one-body stimu- 
lated non-mesonic rates Pi = r„ + Pp for light and medium hypernuclei. 

Other ingredients beyond the OPE might be responsible for the large 
experimental ratios: 

1. Calculations with AA^ — > A^A^ transition potentials including heavy- 
meson-exchange (e.g. K) or direct quark contributions have improved 
the situation. [HI EI 
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2. The analysis of the ratio Tn/Tp is influenced by the two-nucleon in- 
duced process ANN — > NNN: by assuming the quasi-deuteron ap- 
proximation for the absorption of the meson emitted in the A decay, the 
three-body process are mainly Anp nnp and a considerable fraction 
of neutrons could come from this channel in addition to An — > nn and 
Ap np. However the inclusion of the new channel would bring to ex- 
tract from the experiment even larger values for the Tn/Tp ratios. [HI El 

3. The effect of the final state interaction (FSI) on the spectra of the 
emitted nucleons: the nucleon energy /momentum distributions have 
been calculated by using a Monte Carlo simulation to describe the 
nucleon rescattering inside the nucleus. The main effects thus obtained 
indicate that the nuclear collisions remove nucleons from the high to 
the low energy part of the spectrum, moreover the numbers of emitted 
protons and neutrons tend to become similar (due to charge exchange 
processes. 

A recent upgrading of the calculations quoted above shows that the experi- 
mental proton spectra are compatible with values of the T^/Tp ratio between 
0.5 and 1.0, still leaving some discrepancy with the direct theoretical esti- 
mates. 

5 Conclusions 

It is by now well understood that beyond the mesonic channel, hypernuclear 
decay proceeds through non-mesonic processes, induced by one nucleon or 
by a pair of correlated nucleons. This channel is dominant in medium-heavy 
hypernuclei, where the Pauli principle strongly suppresses the mesonic decay. 

The mesonic rates have been reproduced quite well by calculations per- 
formed in different frameworks. The non-mesonic rates have been consid- 
ered within several phenomenological and microscopic models, most of them 
based on the pion exchange. More complex meson exchange potentials and 
direct quark models have also been used for the evaluation of non-mesonic 
decay rates. The obtained rates appear to be in agreement with the exper- 
imental data. 

Although several calculations reproduce the total non-mesonic width, 
Tnm = Tn + rp(+r2), the obtained Tn/Tp are often in strong disagreement 
with the measured central data. Hence further efforts (especially on the 
experimental side) must be invested in order to understand the detailed 
dynamics of the non-mesonic decay. 
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Recent experiments at KEK |19j have considerably reduced the error 
bars on r^/Fp, by means of single nucleon spectra measurements. Good 
statistics coincidence measurements of nn and np emitted pairs are required. 
The angular correlation measurements, as expected from the forthcoming 
FINUDA 20 experiment, will also allow for the identification of nucleons 
coming out from the one- and two-nucleon induced processes. 
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